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We investigate the interaction of Lya photons produced by the first stars in 
■ the universe with intergalactic hydrogen prior to reionization. The background 

Lya spectral profile is obtained by solving a Fokker-Planck equation. Accurate 
O 1 values of the heating and scattering rates, and the spin-kinetic temperature cou- 

pling coefficient, are presented. We show that the heating rate induced by the 
Lya scatterings is much lower than found previously, and is basically negligible. 
The dominant heating source is most likely the X-rays from the first ionizing 
^ ■ sources, which are able to penetrate into the atomic medium. The scattering of 

Lya photons couples the hydrogen spin temperature to the kinetic temperature. 
If the first ionizing sources in the universe did not emit significant X-rays, the 
spin temperature would be rapidly brought down to the very low gas kinetic 
temperature, and a 21cm absorption signal against the CMB larger than 100 mK 
would be predicted. However, we argue that sufficient X-rays are likely to have 
been emitted by the first stellar population, implying that the gas kinetic tem- 
perature should rapidly increase, turning a reduced and brief absorption signal 
into emission, with a smaller amplitude of ~ 10 mK. The detection of the 21cm 
absorption and emission feature would be a hallmark in unravelling the history 
of the "dark age" before reionization. 



Subject headings: cosmology: theory — diffuse radiation — cosmic microwave 
background — intergalactic medium — line: formation — radio lines: general 
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1. Introduction 

After recombination of the primordial plasma at z ~ 1000, the Universe entered an 
age of darkness. The primordial density fluctuations continued to grow, until non-linear 
gravitational collapse led to the formation of the first dark matter halos in which the gas 
could collapse, and then radiatively cool. The dark ages ended at the time when the first 
stars formed in these halos, with a characteristic mass of ~ 1O 6 M (Couchman & Rees 1986; 
Tegmark et al. 1997; Abel et al. 2000; Bromm et al. 2002). These first stars, as well as any 
black holes or neutron stars formed subsequently which could accrete surrounding matter, 
produced photons capable of ionizing hydrogen (with wavelength A < 91.2 nm). As soon as 
the ionizing radiation produced in the sites of star formation was able to escape into the low- 
density intergalactic medium, the process of reionization started. During reionization, the 
ionized volume fraction of the universe increased gradually, until all the low-density regions 
were ionized. Recent observations of the Gunn-Peterson trough (Gunn & Peterson 1965) 
at z > 6 (Becker et al. 2001; Fan et al. 2002, 2003), indicate that the end of reionization 
probably occurred at z ~ 6, although the epoch at which reionization began is still unknown 
(e.g., Gnedin 2000; Miralda-Escude 2003). 

The redshifted 21 cm line opens up a promising window for observing this period of 
cosmic history. This line is produced by the transition between the singlet and triplet 
hyperfine levels of the hydrogen atom at the electronic ground state. The line can be observed 
in absorption or emission against the Cosmic Microwave Background (CMB). The brightness 
temperature is given by 

T s - T C mb n , n noc . -rr \ ( Qbh \ ( 0.3 \ 1/2 (I + z\ 1/2 p HI T s - T C mb 



ST = — — — (1 — e~ T ) ~ (0.025 K) . . , , 

i + z y 1 v ; vo.o3y \n rn0 j v io y p H t s 

(i) 

where T s and Tqmb are the hydrogen spin temperature and the CMB temperature, r is 
the optical depth, pm is the density of atomic hydrogen and pu is the mean hydrogen 
density of the universe (we have used the approximations r <C 1 and f2 m o(l + z) 3 ^> 1). 
Therefore, the detailed spectrum of the CMB at the long wavelengths of the redshifted 
21cm, when observed with high angular resolution, contains precious information on the 
patchwork of density, spin temperature, and ionization inhomogeneities during the era when 
the intergalactic medium was still atomic (Hogan & Rees 1979; Scott & Rees 1990; Kumar, 
Padmanabhan, & Subramanian 1995; Madau, Meiksin, & Rees 1997; Shaver et al 1999; Tozzi 
et al. 2000; Iliev et al. 2002). 

The hydrogen spin temperature plays a key role in determining the amplitude of the 
21cm emission or absorption signal against the CMB. At the same time, the evolution of the 
spin temperature depends on the gas kinetic temperature, and the way that the hyperfine 
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structure level populations couple to the CMB and the kinetic temperatures. The adiabatic 
cooling of the hydrogen due to the expansion of the universe causes the gas temperature to 
decline as (1 + z) 2 , while the CMB temperature is decreasing only as 1 + z, so the gas cools 
below the CMB. At redshifts z > 100, the residual ionization left over from the epoch of 
recombination (Peebles 1968) keeps the gas kinetic temperature at the CMB temperature 
owing to the heating produced by electron scattering of the CMB, but below this redshift 
this heating source becomes small and the gas temperature drops below the CMB. If the spin 
temperature can be coupled to the kinetic temperature at this stage, the atomic medium 
would be observed in absorption. 

The coupling of the spin and kinetic temperatures can be achieved by atomic collisions 
or by scattering of Lya photons. The coupling due to atomic collisions is small except at very 
high redshift, where observations of the redshifted 21cm radiation becomes more difficult, or 
for unrealistically high values of (Scott & Rees 1990). A stronger coupling can be induced 
by Lya photons, which should be produced as soon as the first stars are formed. The Lya 
photons couple the spin and kinetic temperatures by being repeatedly scattered in the gas. 
The resonance scattering of a photon consists of a transition from n — 1 to n — 2, followed 
by the opposite downward transition. This scattering process can change the population 
of the hyperfine structure levels (Wouthuysen 1952; Field 1958). The probabilities for a 
transition back to the ground or excited hyperfine structure states depend on the slope of 
the radiation spectrum near the Lya line center, or the "color temperature". When the 
radiation spectrum reaches a steady state, this color temperature near the Lya line center 
is equal to the gas kinetic temperature, so the Lya scatterings introduce a thermal coupling 
between the spin and kinetic temperature. The result of these couplings is that the spin 
temperature is given by (Field 1958, 1959) 

rp T C mb + y a Tk + y c Tk /0 n 

1+Va + Vc 

where 

Vol = ~, Tf'i He = ~ A 7f ■ \d) 

Here, Tj, is the kinetic temperature, T* = 0.0682K is the hyperfine energy splitting, A w = 
2.87 x 10~ 15 s _1 is the spontaneous emission coefficient of the 21cm line, and C\q = ^iqUh is 
the collisional de-excitation rate of the excited hyperfine level. The value of k\q ranges from 
2 x 10~ 14 to 2.5 x 10 _10 cm 3 s- 1 for T in the range 1 - 1000 K and was tabulated in Allison 
& Dalgarno (1969). 

The indirect de-excitation rate P\o of the hyperfine structure levels is related to the total 
Lya scattering rate P a by P w = 4P a /27 (Field 1958), which is given by P a = f ducn u cr(v), 
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where n u is the number density of photons per unit frequency, and <j{v) is the cross section 
for Lya scattering (Madau, Meiksin, & Rees 1997, hereafter MMR). 

Besides coupling the spin temperature and kinetic temperature of the gas, the scattering 
of Lya photons may also heat up the gas. In this paper, we present the scattering rate and 
heating rate due to Lya photons in a uniform atomic intergalactic medium as a function 
of the two parameters they depend on: the gas temperature and the Lya scattering optical 
depth of the medium (or Gunn- Peterson optical depth). 

There are two different ways in which photons propagating through the atomic medium 
may reach the Lya line, and start being repeatedly scattered by hydrogen: (a) Photons that 
are emitted between the Lya and Ly/3 wavelengths by sources such as massive stars will be 
redshifted as they penetrate the atomic medium until they reach the Lya wavelength. These 
photons enter the Lya line from its blue wing, and will be called continuum photons in this 
paper, (b) Photons that are emitted at wavelengths shorter than Ly/3, but still longer than 
the Lyman limit so that they can penetrate deep into the atomic medium, will redshift until 
they reach the Ly/5 or a higher Lyman series line, and will then be converted into a Lya 
photon when the hydrogen atom they excite decays first to the 2p state before reaching the 
ground state. These photons will be injected near the Lya line center with a distribution 
determined by the Voigt function, and will be called injected photons. These two types of 
photons result in different scattering and heating rates, so results will be presented separately 
for both of them. 

MMR estimated the heating rate due to Lya photons by assuming that the average 
relative change in a Lya photon energy in each scattering is the same as when the photon is 
scattered by an atom at rest: (AE/E) = —hv a /mHC 2 . Hence, they assumed that the total 
energy transfer rate is 

A pi 

E = -(—)hv a P a . (4) 

At the thermalization rate P t h = 27 Ai Tcmb/^T^ (which is the rate P a required to bring 
the spin temperature close to the kinetic temperature), the heating rate assumed by MMR 
is E/kB — 220K Gyr _1 (l + z) /7, which would heat up the gas kinetic temperature above the 
CMB temperature in a fraction of a Hubble time at z ~ 6. 

In the present investigation we calculate the heating rate of the continuum and in- 
jected Lya photons by computing the background spectrum near the Lya line using the 
Fokker-Planck method introduced by Rybicki & Dell'antonio (1994) (hereafter RD), and 
using energy conservation. We find that the heating rate due to Lya scattering is much 
smaller than the estimate of eq. (4), and as a consequence Lya photons do not constitute a 
significant heating source of the gas. Calculating the heating rate caused by Lya photons as 
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if every atom were at rest, without taking into account the thermal motions, is incorrect, and 
as we shall see the heating rate per unit volume is not proportional to the atomic density. 
We explain in the Appendix the reason why the argument used by MMR in their Appendix 
B to justify using eq. (4) for the heating rate is invalid. As pointed out by MMR, X-rays are 
important in pre- heating the atomic medium before it is reionized. In fact, our work leads 
to the conclusion that X-rays are the only important source of heating before the atomic 
medium is reached by an ionization front. 

This paper is organized as follows: in §2 we describe our method of calculation. We 
present our result on scattering and heating rates in §3, explaining the physical reason why 
gas is not heated by Lya photons as if every atom were at rest. We illustrate the consequences 
for future 21cm observations in §4, with an example for the evolution of the spin temperature 
with redshift, including also a model for X-ray heating. Our conclusions are summarized in 
§5. We adopt a flat ACDM model with the following cosmological parameters: h = 0.7, 
tt A = 0.7, Vt m = 0.3, tt b h 2 = 0.021, Y Hc = 0.24. 

The code for solving the line profile and calculating the scattering and heating rates is 
available from the authors 1 . 



Our goal in this work is to compute the Lya scattering rate and the heating rate of 
hydrogen induced by this scattering, from which one can then compute the evolution of the 
spin temperature. The scattering and heating rates can be calculated once the background 
spectrum near the Lya line is known. Let us consider the background of photons produced 
by constant, flat spectrum sources in an expanding universe (in practice the intrinsic source 
spectrum does not matter because Lya scattering is important only on a very narrow range 
of frequencies around Lya , in which the intrinsic source spectrum can be considered to be 
constant). In the absence of scattering, the photons would simply redshift and would produce 
a photon number density per unit frequency n(u) = no that is independent of frequency near 
Lya . The presence of scatterings by hydrogen introduces a feature in the spectrum around 
the Lya line, which has a constant shape n{v) once a steady-state is reached. Given this 
spectrum, the rate at which Lya photons are being scattered is given by the integral 



1 Our code, named "Lyman Alpha scattering and Spin Temperature", or LAST, can be downloaded at 
http://theory.kitp.ucsb.edu/~xuelei/LAST. 



2. 



Method of Calculation 




(5) 



- 6- 



Here, cr(y) is the Lya scattering cross section as a function of frequency. The cross section 
has a thermal width Av D determined by the temperature T of hydrogen atoms of mass m# , 



^2k B T/m H 

Av D = v a , (6) 

c 

where z/ Q is the Lya line central frequency and ks is the Boltzmann constant. The cross 
section at the line center is 

(7 = /^[A^]- 1 , (7) 

m e c 

where /12 = 0.416 is the oscillator strength for the Lja transition. The dependence of the 
cross section on frequency is given by the normalized Voigt function, 

a f°° e~ y2 
* {x) = l?&L dy {x-yy + a" (8) 

where x = (v — v a )/AvD, the Voigt parameter is a = A 2 \/ {SnAun), and A 2 \ = 6.25 x 10 s 
Hz is the Einstein spontaneous emission coefficient. We will also use the photon intensity 
J(u) = cn(v)/(4ir). With these definitions, the scattering rate is 

P a = 4ira Av D J J(x)(f>(x)dx . (9) 

The heating rate of the gas induced by Lya scattering can be simply obtained by 
requiring conservation of energy. The scatterings cause the cosmic background to have a 
narrow spectral feature around the Lya wavelength. Were the scatterings to cease at a certain 
instant, this feature would simply be redshifted to longer wavelengths. Instead, when a 
steady-state spectrum is reached, the feature is fixed at a constant wavelength. The reduction 
in energy per unit volume of the background due to this feature is Ae = J dv[riQ — n(u)] hv dv 
(for the purpose of this discussion we ignore the photons injected at the Lya wavelength by 
recombinations, so that in the absence of scatterings, n{y) = no; injected photons will 
be fully considered later). Over an interval of time dt, the energy per unit volume that 
the background photons must lose to keep the spectral feature at a fixed wavelength is 
Ae{dv/v) = AeH dt, where dv is the change in frequency due to redshift over the interval 
dt, and H is the Hubble constant at time t. The energy lost by the photons must be gained 
by the hydrogen atoms, so the gas heating rate per unit volume is 

T = HAe = H J [no- n{y)\ hv dv = 4nHh ^ AuD J [j Q _ j^)] v dv . ( 10 ) 

We use these equations in §3 to compute the scattering and heating rates once the line 
spectrum has been calculated. An alternative derivation of the heating rate in equation (10) 
is given in the Appendix. 
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The spectrum of the cosmic background near Lya is determined by the radiative transfer 
equation for resonant scattering for a homogeneous and isotropic expanding universe, which 
is (RD) 

^^ = -^ x )J^r)+ ls ^^- + J Rn(x, x') J(x', r)dx' + C{r)^{x) , (11) 

where r = (cnicr )t is the mean free time between scatterings for a photon at the line center, 
t is the physical time, and ri\ ~ n# is the number density of hydrogen atoms in the ground 
state. The Sobolev parameter is 75 = rjp, where tqp = (11 e 2 ri\ f 12) / (Hm e v a ) is the Gunn- 
Peterson optical depth. The first term in the right-hand-side of equation (11) describes 
the removal of photons that are scattered, and the second term accounts for the redshift. 
The third term represents the addition of the photons redistributed over frequency after 
scattering (Rybicki & Hummer 1992), and the fourth term is due to the injection of new 
Lya photons with a spectral profile ip(x). 

In the Fokker-Planck approximation, the scattering problem is treated as diffusion in 
the approximation that the width of the spectral feature is large compared to the frequency 
change in typical scatterings. The redistribution term then becomes (RD) 



/ 



1 d 

R H (x,x')J(x' \t)dx' = Mx)J(x) + - — 

2 ox 



^)^+2#(x)J(x) 



(12) 



where the recoil term is included, with rj = hv a j (c\ / 2rriHkBT). This approximation conserves 
photon number. The equation is then reduced to [see eq. (34) of RD] 



dJ_ _ ld_ 

dr 2 dx 



dJ 



(97 

+ -Ys-fo+C(t)1>(x). (13) 



When the equilibrium state is reached, dJ/dr = 0. 

Since this equation is linear, we can write its solution as J = J c + Jj, where J c is the 
spectrum of the continuum photons, and Jj is the spectrum of the recombination photons 
injected near the line center. 

For the continuum photons, there is no injection and so C(i) = 0. Equation (13) is 
reduced to 

<j)(x) — J c (x) + 2 [ V (f>(x) + 75] Jc(x) = A, (14) 

where A is a constant of integration determined by the boundary condition: as x — > 00, 
4>(x) — > 0, J c — > J c o, where J c o = (47r/c)n c o is the unperturbed intensity. Therefore, 
A = 27 5 J . The equation is then solved as an initial value problem of an ordinary dif- 
ferential equation with the condition J c (x) — > J c0 for x — > —00, required by photon number 
conservation. 
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In Fig. 1, we plot the solution for a typical case: T = 2.65K and 75 = 1.44 x 1CT 6 , 
corresponding to the physical conditions of unheated gas at z = 10 in a flat ACDM model. 
An absorption feature is produced near the center of the line. 

For the injected photons, we assume a constant source C(t) = C. This should generally 
be an accurate approximation because the timescale for the evolution of the production rate 
of Lya photons from sources is likely be of order the Hubble time, which is much longer 
than the relaxation time, which is of order the time to redshift across the width of the Lya 
spectral profile of the cosmic background. We also inject photons with a single frequency at 
the Lya line center, i.e., ip(x) = 5(x). Because photons change randomly in frequency over 
a thermal width at each scattering, and the number of scatterings each photon undergoes is 
very large, injecting the photons with a thermal width instead does not make any difference 
to the final profile [we have verified this by obtaining solutions with tp(x) given by a Voigt 
profile, which did not result in any significant change to the resulting spectral profile]. The 
solution is 

(f>(x) -^Hx) + 2 [#(a;) + 7s] Ji(x) = -2CG(x) + B , (15) 

where Q(x) is the Heaviside step function. The value of the integration constant B can again 
be determined from boundary conditions: as x — * +00, we must now have Ji(x) — > 0, from 
which we obtain B = 2C. The profile of the line can again be obtained by solving eq. (15), 
with the initial condition Jj(x = —00) = J«o, where J i0 is the photon intensity due to the 
injected photons. The result is shown in Fig. 2, for two different temperatures; we see that 
some of the photons injected at the line center are scattered to the right (blue) side. 

We note that the solution for the spectral profile of continuum and injected photons 
has a characteristic width in the variable x = u/Aud that is substantially larger than unity, 
implying that the Fokker-Planck approximation we have used should be accurate. 



3. Results 

The heating and scattering rates are of course proportional to the background intensity. 
During the epoch of reionization, at least one ionizing photon has to be emitted for every 
baryon in the universe, and more if the typical baryon can recombine many times. The 
number of Lyct photons emitted is of the same order for sources without a large Lyman 
break, and can be much larger if many ionizing photons are internally absorbed. We use the 
intensity corresponding to one photon per hydrogen atom in the universe as a fiducial value: 

Jo(*) = ^ = 1-65 x 10- 13 (1 + zf (^f) cm - V i sr- Hz" 1 , (16) 
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where nn is the hydrogen number density. For reference to other work, we quote also the 
value ofthe energy intensity of this unit: Johu a = 2.69xl0~ 24 (l+z) 3 (£l b hl/0.02) erg cm -2 s _1 sr -1 Hz" 
The heating rate in equation (10) can be reexpressed as a heating rate per hydrogen atom 
and per Hubble time. For the continuum photons, it is: 

T c hAu D J c0 



Hn H k B k B J 



[dx l- J -M. ^^n% J -Sl C) (17) 

J \_ -JcO J Jq 



where we have defined I c = J dx[l — J c (x)/J c q], and T = {hv a ) 2 /(k B mHC 2 ) = 1.29 x 10~ 3 
K. For the example shown in Figure 1, the integral is the shaded area, I c = 16.7. It is clear 
from equation (17) that if the number of photons in the background is not greater than one 
per hydrogen atom, the heating rate due to Lya scattering is not very large because To <C T. 
In the example of T = 2.65 K at z = 10, equation (17) yields T c / '(Hnnks) = 1.4 K (J c0 / 'Jo). 

Similarly, for the injected photons with a 5-function profile, the heating rate is 

Hn H k B k B J [J-oo L Jio J Jo Jio J Jo 

We note that, for the more exact case in which photons are injected with a Voigt profile in 
frequency, the heating rate is obtained by subtracting from Ji(x) the steady-state solution 
of the background with no scattering (shown in Fig. 3 of RD). However, one can easily 
prove that the result is exactly the same as what is obtained by subtracting instead the step 
function spectrum that results from a 5 function of the injected photons. 

In Fig. 2, we show the spectrum Jj for two different temperatures: = 2.65K and T = 
IK, for 75 = 1.44 x 10~ 6 . In the T k = 2.65K case, we have Ti/(Hn H k B ) = -0.104K(J i0 /J ). 
This negative heating rate means that the scattering of these injected photons actually cools 
the gas, because some of the photons are scattered to the blue side of the spectrum and gain 
energy from the hydrogen gas. At lower gas temperatures, as in the = IK case shown 
in the figure, Ti/(Hn H k B ) = 0.40 K ( J i0 / J ), the photons heat the gas. However, at z — 10 
the temperature of adiabatically expanding gas is 2.65 K, and there is no known mechanism 
which could cool the gas below this temperature, so the net effect of the injected photons 
must be to cool the gas. In any case, both the heating rate due to continuum photons and 
the cooling rate due to injected photons are very small. 

To calculate the spin temperature from equation (2), we also need to know the scattering 
rate of the Lya photons in equation (9), which can be reexpressed as 

P a = Hr GP [ ^-^dx. (19) 
J Jo 



We define 



S c = f ^^(x) dx, Si= [ J -^4>{x) dx. (20) 

J JcO J JiO 
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We plot the integrals I c , I i: S c , and Si as a function of 75 = r^p, for several different 
temperatures, in Figures 3- 4. Our code calculates the values of these integrals for any given 
gas temperature and optical depth. 

The intensity corresponding to the thermalization rate required to bring the spin tem- 
perature down to the kinetic temperature (see the introduction) is 

^0 _ Pa _ 27A w T C mb1s /«-,-, x 

Jo ~ Hr GP S A%HS 1 ' 

where the subscripts c or % can be applied. For the case 75 = 1.44 x 10 -6 and T = 2.65 K 
at z — 10, Si or S c are close to one and we find Jo/ Jo ~ 0.3. Thus, from our earlier result 
that the heating rate in this case is T c /(Hn H k B ) = 1.4 K (J c o/Jo) for the case of continuum 
photons, we find that over the time that continuum Lya photons bring the spin temperature 
down to the gas temperature, they should heat the gas by only ~ 0.5 K, and injected photons 
should cool the gas by ~ 0.03 K. 

The reason why the heating rate is much less than the value obtained by assuming 
that each atom loses the same average energy at each scattering as if it were at rest is that 
the intensity of the background spectrum drops steeply at the Lyct line center, as shown in 
Figures 1 and 2. The vast majority of the scatterings occur within a few thermal widths of 
the line center, where the color temperature of the spectrum (reflecting its slope) is almost 
equal to the kinetic temperature. When these two temperatures are the same, there obviously 
cannot be any heat exchange between photons and atoms because they are in thermodynamic 
equilibrium. The heating from continuum photons occurs because the color temperature is 
slightly larger than the kinetic temperature. Injected photons have a temperature slightly 
lower than the kinetic temperature, so they cool the gas. Cooling is caused by photons that 
are on the red side of the line center because these photons are more likely to be scattered 
by atoms moving in the direction opposite to the photon (for which the Doppler effect can 
shift the photon frequency to the line center), and these atoms will be slowed down by the 
scattering. When the color temperature is below the kinetic temperature, the cooling effect 
due to the preferential scattering of photons moving in the opposite direction than the atoms 
is greater than the heating effect due to recoil from photons moving in a random direction. 

The heating or cooling rates per unit volume are not proportional to the atom density. 
If the atom density is increased, not only the number of scatterings increases, but also the 
color temperature is brought closer to the kinetic temperature because of the increased rate of 
interaction between atoms and photons, and the average energy exchanged between them at 
each scattering is reduced. The heating rate changes only with the width of the background 
spectral feature that determines the integrals I c and according to the energy conservation 
argument used in §2. 
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In the calculation of the spin temperature evolution from equation (2), we assume that 
the color temperature of the Lja photons is equal to the gas kinetic temperature. This is 
true at the line center in the limit of large optical depth (Field 1958, 1959), but the slope of 
the photon distribution varies outside the line center. In the more general case, equations 
(2)- (3) can be generalized as 

y a T k -> Jdx y a (x)T a (x) = J^^J dP a (x) (22) 

Va^ Jdx y a (x) = J dP a {x)/T a {x) (23) 

where the color temperature T a at x is defined by 

J'{x) _ hv hu 



(24) 



J{x) k B T a (x) k B T a (x) ' 
and 

dP a (x) = 4ira Au D J(x)(f)(x)dx (25) 

Because the vast majority of scatterings occur near the line center, where the color tem- 
perature and kinetic temperature are very nearly the same, this generalization does not 
significantly change the result for the spin temperature. We have found by a numerical test 
that this correction affects the spin temperature by only a small fraction of a percent. 



4. Simple Models for the Thermal Evolution 

We will now consider a simple model where the emissivity of Lya photons turns on at a 
redshift Zi and increases linearly with redshift thereafter. We will also consider the effect of 
X-ray emission from the same star formation regions, which dominates the heating rate of 
the gas. This will illustrate a plausible thermal history of the gas and the expected strength 
of the absorption or emission signal in the redshifted 21cm line. 

The kinetic temperature of the gas in the expanding Universe evolves as 

h -u -y\ dTk - 9T 2 (T tot - A tot ) 

(1 + ^- 2Tfe ~3 (Jinks) ' (26) 
where T tot and A tot are the total heating and cooling rates, respectively, and n is the total 
gas number density. 

In the absence of heating and cooling, the gas temperature decreases adiabatically with 
T oc (1 + zf. We use the code RECFAST 2 (Seager, Sasselov, & Scott 1999; Seager, 



2 http: / / www-cfa.harvard.edu/~sasselov/rec 
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Sasselov & Scott 2000) to calculate the temperature evolution of the gas before the first 
stars and quasars turn on. We note that at z > 40 the spin temperature drops below the 
CMB temperature due to the collisional coupling (in Fig. 6, this drop is seen at the highest 
redshift end). The corresponding absorption is, however, difficult to observe at the very 
long wavelengths corresponding to this redshift. At lower redshifts, the spin temperature is 
practically equal to the CMB temperature until the first Lya photons are produced. 



4.1. Lya emissivity 

If the comoving photon emissivity (defined as the number of photons emitted per unit 
comoving volume, time and frequency) at the Lya line from stars at redshift z is e(u a ,z), 
then the comoving Lya photon intensity at redshift z is given by 



■maa: 



■Uz) = 1(1 + zr jf H^Tf^ + z')- (1+as) dz', (27) 

where the emissivity spectrum is assumed to be e(v) oc z/~ as_1 , and z max is determined by 
the condition that only the photons emitted at frequencies lower than the Ly/5 line will not 
interact with the intergalactic medium: 

(1 + z max ) = ^(1 + z) = — (1 + z) (28) 
v a It 

To discuss the plausible range of values of the comoving emissivity, it is convenient to 
introduce the unit 

i=?** = *!LHJ . (29) 

Vol C 

Thus, if e = e, one photon is being emitted per hydrogen atom and per Hubble time, per 
unit \og(v). 

During the epoch of reionization, at least one photon must be emitted for every baryon 
in the universe. For a homogeneous medium, the recombination rate om\ is equal to the 
Hubble rate, H, at z ~ 6; hence, if reionization takes place during the interval 20 > z > 6, 
a few photons should be emitted for every baryon to ionize them and to compensate for 
recombination. The number of recombinations can be increased due to the clumping factor 
of the ionized intergalactic gas (see, however, Miralda-Escude, Haehnelt & Rees 2000 for a 
discussion of why the clumping factor of ionized gas is probably not much larger than unity 
during the epoch of reionization). If these ionizing photons are produced by star-forming 
regions, then the photon emissivity per unit log(z/) between Lya and the Lyman limit should 
be 3 to 5 times larger than the emissivity of ionizing photons, owing to the average Lyman 
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break of the atmospheres of massive stars with a Salpeter Initial Mass Function (e.g., Bruzual 
& Chariot 1993). The photons between Lya and the Lyman limit may have a larger escape 
fraction from the star formation regions than the ionizing photons (if some of the absorption 
is due to local hydrogen with low dust content), in which case the ratio of the emissivity 
between Lya and Lyman limit to the ionizing emissivity would be further increased. This 
suggests that the emissivity near Lya should be at least e > lOe. 

The total number of ionizing photons emitted is also related to the metallicity produced 
by massive stars. To produce one ionizing photon per baryon, the mean metallicity of the 
universe must be raised by 10 -5 , or 5 x 10 _4 Z Q (e.g., Madau & Shull 1996). If the typical 
metallicity in the Lya forest and damped Lya systems at z > 3, which is Z ~ 10~ 2 Z Q , is not 
to be exceeded, one cannot produce more than ~ 20 ionizing photons per baryon, implying 
that e < lOOe at Lya . 

As a simple model of the Lya background evolution, we assume that the Lya emissivity 
increases linearly from at some initial redshift Zi, 

i \ I 0j z > z i (m} 

~ I € rfE^\ Z r <Z< Zi 

We will use two values for the peak emissivity at Lya , e r , reached at z r : e r = lOe and 
e r = lOOe, based on the previous discussion. Assuming also a flat spectrum, as = 0, the 
Lya intensity due to the continuum photons is then 

-g{i + z) , (31) 



J e 



where the function g is 



x 3/2 

g( X ) = — 

Zi z r 

and 



1 + (V 3 / 2 - x- 3/2 ) - 2 (V 1 / 2 - x- 1/2 ) 



x<- = mm 



z<Zi (32) 



(33) 



The resulting proper Lya intensity is plotted as the function (1 + z) 3 J c0 / J in Fig. 5, for 
z r = 6 and e r /e = 10, for two cases of initial redshift: Zi = 20 and Zi = 12. The proper 
Lya intensity increases rapidly after Zi and then decreases as z decreases due to the (1 + z) 3 
factor. 



We will assume that J i0 = J c0 . The injected photons have a comparable intensity to the 
continuum ones, because the photons emitted between Ly/3 and the Lyman limit are similar 
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to those emitted between Lya and Ly/3, and a large fraction of the Lyman continuum photons 
will eventually result in an injected Lya photon when the atom that they ionize recombines 
again. 



Since the scattering of Lya photons produces a negligible heating rate for the gas, the 
main heating mechanism should be soft X-ray photons. The soft X-rays can be emitted by 
a number of sources associated with star-forming regions, such as X-ray binaries, supernova 
remnants and QSOs. 

X-ray photons have a photoionization cross section that is much smaller than that of the 
ionizing photons near the hydrogen ionization potential, so they can penetrate deep into the 
neutral regions of the intergalactic medium. There, they photoionize hydrogen or helium, 
and the high-energy electron that is produced can then transfer a fraction of its energy to 
the gas by Coulomb collisions with other electrons and protons. Once the energy has been 
thermalized among the charged particles at a temperature low enough to avoid losses by 
collisional excitation of hydrogen, the thermal energy can be transferred to the hydrogen 
atoms in slow collisions involving dipole interaction. The fraction of the X-ray energy that 
can be converted to the gas thermal energy (rather than being used in collisional ionizations 
and excitations of hydrogen or helium) depends on the fractional ionization of the gas. Using 
the RECFAST code, we found that the initial ionized fraction from the recombination era 
(before more ionizations are caused by the X-rays) is 2 x 10 -4 , which implies that a fraction 
of 14% of the X-ray energy is used to heat the gas (Shull & Van Steenberg 1985). As the 
fraction of free electrons increases, this fraction also increases, but for simplicity we shall 
assume it a constant. 

We parameterize the emissivity in X-rays in terms of the fraction of the energy that is 
emitted in X-rays compared to the energy emitted at Lya per unit log u, which we designate 
as a x . The X-ray heating rate is then 



where e(z) is given in Eq. (30), and h is the Planck constant. At z — 6, the heating rate due 
to X-rays in this model is 



Therefore, by emitting only ~ 1% of the energy as X-rays, the sources can rapidly heat the 
intergalactic medium to temperatures above the CMB. 



4.2. Heating by soft X-rays 




(34) 




(35) 
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4.3. Results for thermal evolution and the 21 cm Antenna Temperature 

The evolution of the gas kinetic and spin temperatures in the absence of X-ray heating 
is plotted in Fig. 6, for two values of the initial redshift (zi = 20 and Zi = 12), and two values 
of the emissivity [e r je = 10 and 100). 

In the absence of X-rays, the heating due to Lja scattering is very low and the kinetic 
temperature of the gas deviates only slightly from the adiabatic curve at late times. For our 
assumed emissivity of 10 to 100 photons per baryon per Hubble time, the spin temperature 
drops to values very close to the kinetic temperature, producing a large absorption signal on 
the CMB of up to 200 mK, as seen in Figure 7. 

Figure 8 shows the thermal evolution and the resulting antenna temperature when we 
include X-ray heating with a x = 0.01, for the case Zi = 12 and e r /e = 10. With only 1% 
of the energy emitted as X-rays, the kinetic temperature of the gas increases rapidly when 
the first sources turn on, leaving only a brief period of absorption in the beginning of small 
amplitude, which then turns to emission, with an amplitude of less than 20 mK. 

5. Conclusion 

We have calculated the Lja photon scattering and heating rates in the atomic inter- 
galactic medium, including both redshifted continuum photons emitted between Lja and 
Ly/3, and injected photons produced by recombinations or excitations by Ly/3 and higher 
Lyman series photons. The heating rate of the gas due to the scattering of these Lja pho- 
tons is very small and essentially negligible. In the absence of other heating sources, the gas 
would continue to cool down almost adiabatically while the spin temperature would rapidly 
drop to nearly the kinetic temperature given a reasonable emission of photons during the 
reionization epoch. This would produce a 21cm absorption against the CMB as large as ~ 
200 mK, at the average hydrogen density of the universe. 

Of course, if the absorption signal can be observed at high angular and frequency res- 
olution, the antenna temperature should fluctuate around its mean value owing to density 
fluctuations, as well as spin temperature fluctuations that would be produced by spatial vari- 
ations in the history of the photon emissivity between Lja and the Lyman limit (e.g.Tozzi 
et al. 2000; Ciardi & Madau 2003). At the same time, during the epoch of reionization, there 
should be a gradual increase of the fraction of the volume that is ionized (e.g., Miralda-Escude 
2003), and the average antenna temperature should then be multiplied by the fraction of 
baryons that remain in atomic form. Because this fraction of baryons should decrease grad- 
ually, any 21cm absorption or emission should also vary smoothly with redshift. 
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Unfortunately, the large absorption signal predicted in the absence of X-ray heating is 
unlikely to be present. It is very likely that significant X-ray emission is produced within 
only a few million years of the formation of the first stars in the universe: some of the stars 
may be in binaries which will evolve into high-mass X-ray binaries, and supernova explosions 
can also emit soft X-rays (see e.g. Oh 2001; Venkatesan, Giroux, & Shull 2001). Photons 
with energies as low as 0.2 keV can propagate as far as 3 comoving Mpc over the atomic 
intergalactic medium at z ~ 10 before they are absorbed by hydrogen or helium. 

There may still be some ways in which certain regions of the intergalactic medium may 
be illuminated by Lya photons from the first stars, without receiving soft X-rays as well. 
For example, if stars form in a region with a high column density of hydrogen and helium 
with very little dust, and the hydrogen and/or helium is not fully ionized for a long time, the 
photons between Lya and the Lyman limit from stars would escape, while the ionizing and 
soft X-ray photons might all be locally absorbed. In particular, there is the possibility that 
singly ionized helium remains self-shielded as massive stars ionize the gas around them. We 
note from Figure 6 that the time required to bring down the spin temperature can be very 
short (10 8 years or less) for typical values of the Lya background intensity, and it is possible 
that a complete absorption of soft X-rays from the first star-forming regions would persist 
long enough to allow for a strong 21cm absorption signal in local regions of the intergalactic 
medium, when they are illuminated by the first light from stars forming in their vicinity. 
Nevertheless, it seems likely that any decrease of the spin temperature will be short-lived, 
and that X-ray heating will cause the 21cm signal to switch to emission a short time after 
the first stars are formed. 

Future observations of the 21cm signal in absorption and emission at high angular and 
frequency resolution have an enormous potential for revealing the details of the reioniza- 
tion epoch. Fluctuations in the antenna temperature arise from both gas density and spin 
temperature variations. The density fluctuations will inform us of the state of gravitational 
evolution of primordial fluctuations at the highest attainable redshifts after recombination, 
while the spin temperature fluctuations will reveal the spatial distribution of the first sources 
of ultraviolet and X-ray emission in the universe. 

This work was supported in part by NSF grants NSF-0098515 and PHY99-07949. 

Appendix 

MMR assumed that the heating rate per unit volume of the hydrogen gas due to Lya 
scattering is equal to the recoil energy given to a hydrogen atom at rest, {hu a ) 2 /(uinc 2 ), 
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times the scattering rate: 

Fmmr = n H P a = kirn H cr Au D / J(x)(f>(x) dx . (1) 

m H c z m H c z J 

This heating rate was justified from equation (13) for the time evolution of the photon 
intensity J. MMR assumed that the rate of heat transfer to the atoms is determined by the 
change in the spectrum due to the recoil effect, which is the second term in the right-hand- 
side of equation (13): 

dJ\ _ djr^J) 

dT J recoil 9X 

They then computed the heating rate by integrating the variation in photon energy density 
per unit physical time over frequency: 

?mmr = ^hAu D n H a J (Ik^ = (Hn H k B )T T GP ^j- S c , (3) 

where the integral S c is given in equation (20), and tqp is the Gunn-Peterson optical depth 
(all symbols are defined in §2 and 3; the last equality involves an integration by parts). 

However, the first term in equation (13) must also be included to compute the heating 
rate of the gas. In fact, the heating rate of the gas is not just determined by the recoil, but 
also by the fact that the background intensity is greater on the red wing of the Lya line 
than on the blue wing, and scatterings have an average tendency to bring photons back into 
the line center. The energy change of the photons due to the first term in equation (13) 
can only go into the gas (whereas that of the third term is due to redshift and goes into the 
expansion of the universe, and that of the fourth term is provided by a source). Thus, the 
correct heating rate is 

dxX dx \ 2 ~dx +r] ^ J ) ' ( 4 ) 

Under steady-state conditions (and no source term; here, we consider continuum pho- 
tons only, to simplify the discussion), we can replace the quantity inside parenthesis by 
—TQpdJ/dx, and integrating by parts, we obtain 

T = (Hn H k B )Vnf^I c , (5) 

which is the expression in equation (17). Note that, because J = (n#c)/(47rz/ Q ), the heating 
rate per unit volume is independent of the atomic density, except for the weak dependence 
of I c on tgp, whereas MMR assumed that the heating rate per atom is independent of n#. 
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X 



Fig. 1. — Continuum Lya photon spectrum. Solid line marks the flat (no-scattering) limit, 
dash line for spectrum produced with T=2.65 K, 75 = 1.44 x 10~ 6 , corresponding to the 
condition of unheated gas at z — 10 in ACDM model. The area of the shaded region gives 
the heating rate. 
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Fig. 2. — Injected Lya photon scattering spectrum. Solid line shows the no-scattering limit. 
The dotted line for Tk = 2.65 K gas at z — 10 in which photon gains energy, dashed line for 
T k = 1 K gas at z = 10 in which photon looses energy. 





Fig. 4. — Scattering rate integral S c (left) and Si (right) versus log 75 for different tempera- 
tures. 
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Fig. 5. — The Lyman alpha proper intensity, plotted as (1 + z) 3 (J c0 /J ), as a function of 
redshift. Dashed and solid curves are for Zi = 20 and Zi = 12, respectively. We also assume 

JiO = JcO 
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Fig. 6. — The kinetic and spin temperature evolution of the gas with only Lja photons 
which increases linearly after some initial redshift Z{. The CMB temperature and adiabatic 
evolution of gas kinetic temperature is plotted in dotted lines, Lya heated kinetic tempera- 
ture in dashed lines, and spin temperature in solid lines. The two curves are for fudge factor 
jfe = 1 and 10. Left: The "late start" scenario with z { = 12. Right: The "early start" 
scenario with Zi = 20. Results for two Lya emissivity models are plotted in each figure, with 
fudge factor ^ = 1 and 10. 




Fig. 7. — The antenna temperature difference as a function of redshift with only Lya photons. 
Left: late start Zi = 12; Right: early start z,- L = 20. 
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Fig. 8. — Model with X-ray heating, Zi — 12, ffe — 1, a x = 0.01. Left: evolution of 
temperatures; Right: the antenna temperature difference. 



